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The phenomenon of synchronization occurring in a locally coupled map lattice subject to an 
external drive is compared to the synchronization process in an autonomous coupled map system 
with similar local couplings plus a global interaction. It is shown that chaotic synchronized states 
in both systems are equivalent, but the collective states arising after the chaotic synchronized state 
becomes unstable can be different in these two systems. It is found that the external drive induces 
chaotic synchronization as well as synchronization of unstable periodic orbits of the local dynamics in 
the driven lattice. On the other hand, the addition of a global interaction in the autonomous system 
allows for chaotic synchronization that is not possible in a large coupled map system possessing only 
local couplings. 

PACS numbers: 05.45.-a, 89.75.Kd 
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Coupled map lattices have provided fruitful theoreti- 
cal and computationally efficient models for investigating 
a variety of processes in spatially distributed dynami- 
cal systems, such as synchronization, pattern formation, 
phase separation, turbulence, nontrivial collective behav- 
ior, etc. . Recently, the phenomena of synchroniza- 
tion and pattern formation induced by external forcing 
on spatiotemporal dynamical systems, such as chemical 
reactions 0- S have been the focus of much attention. 
There has also been interest in experimental investiga- 
tions of spontaneous pattern formation and emergence of 
synchronization in spatially extended systems of interact- 
ing dynamical elements, such as one-dimensional arrays 
of electrochemical oscillators [Bj], chemical and hydrody- 
namical systems with global coupling [|| Q , and popu- 
lations of chaotic electrochemical cells having both lo- 
cal and global interactions The relationship between 
forced spatiotemporal systems and autonomous dynam- 
ical systems possessing global interactions has recently 
been explored in the framework of coupled map lattices 
0. In this paper, we investigate the emergence of syn- 
chronization in forced spatiotemporal systems by using a 
model of a coupled chaotic map lattice subjected to an 
external drive. We show that the chaotic synchronized 
state in this lattice is analogous to the chaotic synchro- 
nized state emerging in an autonomous coupled map sys- 
tem having similar local couplings plus a global interac- 
tion, but other collective states occur differently in these 
two systems. We also show that the addition of a global 
interaction in the autonomous system allows for chaotic 
synchronization that is not possible in a large coupled 
map system possessing only local couplings. 

As a model of a driven spatiotemporal system, we con- 
sider a one-dimensional coupled map lattice subjected to 
a uniform external drive 0, 



x\ +l =(l-e 2 )/(a4) 

+ f [{f{x\ +1 ) + fix*' 1 ) 2f(xt)] + e 2 F t . 



(1) 



scribes the local chaotic dynamics; ei measures the local 
diffusive coupling, e 2 expresses the coupling to the exter- 
nal forcing, and F t is the uniform driving term which can 
be any function of time. 

The dynamics of the driven lattice can be compared 
with that of an autonomous spatiotemporal system, de- 
scribed by the following coupled map system possessing 
both, local and global interactions, 



x\ +l = (1 - e 2 )f{x\) + f + /K" 1 ) 

-2/(*D] + #EL/04), 



(2) 



where f(x) is the same local map as in Eq. ei and 
e 2 are the local and global coupling parameters, respec- 
tively; and the global interaction is provided by the mean 
field of the system. Periodic boundary conditions are as- 
sumed in both systems, Eq. Q and Eq. J5J • 

When the autonomous coupled map system in Eq. (0) 
reaches a synchronized state at some parameter values, 
the evolution of its mean field is identical to the dynamics 
of any local map. Thus for the same set of parameters, 
the driven lattice, Eq. (Q), subjected to a forcing F t equal 
to the local map should exhibit a synchronized state sim- 
ilar to that of the associated autonomous coupled map 
system Eq. © . 

The driven lattice, Eq. can be expressed in vector 
form as 



x t+ i 



(l-e 2 )7+^L]f(x t ) + e 2 F i! (3) 



where x\ is the state of element i 



(i = 1,2,..., N) at 
discrete time t; N is the number of oscillators; f(x) de- 



where the vector components are [x^L = x\, [f(x t )] i = 
f{x\) and [F t ]i = F t ; I is the N x N identity matrix, and 
L is the N x N matrix expressing the diffusive coupling 
among the elements, with components given by La±i = 
1, La — —2, and — 0, otherwise. 

The driven lattice may reach different asymptotic spa- 
tiotemporal patterns, ordered, synchronized, or turbu- 
lent, depending on the characteristics of the drive func- 
tion F t , and on the initial conditions Q . Here we consider 
synchronized states of the lattice induced by a periodic or 
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a chaotic drive. A synchronized state at time t is denned 
by the condition x\ — x tl Vi. In a synchronized state, the 
driven system must satisfy 



x t+1 = (1 - e 2 ) f(xt) + eiF t . 



(4) 



The linear stability analysis of the synchronized state 
in the driven system leads to the condition 



1 



62 + 2 Vj 



< 1 



(5) 



where Vj = -4 sin 2 (nj/N) , j = 0, 1, 2, . . . , N - 1, is the 
set of eigenvalues of the coupling matrix L, with N/2 dis- 
tinct eigenvalues and each being doubly degenerate [lflj : 

and A = limx^oo A X)t=o 1 1°§ l/'C 2 -*)! ^ s the Lyapunov ex- 
ponent of the local map. Thus the range of the parameter 
e 2 for which a synchronized state is stable corresponds to 



l-2e lS in 2 (^)- e - A < e 2 < l-2ei sin 2 (^) + e - A . 



(6) 



In particular, in a synchronized periodic state all the 
elements follow the same cyclic sequence of values. Con- 
sider, for example, an orbit of period p of the local map, 
defined by f^ p '(x n ) = x n , where {x\, a; 2 , ■ ■ ■ , x p } are the 
set of consecutive points belonging to the orbit, satis- 
fying f(x n ) = x n+ i, f(x p ) = xi. This periodic orbit 
is unstable if e A = []^ =1 \f'(x n )\ > 1. If an unsta- 
ble periodic orbit gets synchronized in the driven lat- 
tice, then xt+i — x n +i and xt — x n , and Eqs. (@J yield 
the solution F t = Thus, if the external drive 

follows a periodic unstable orbit of the local map, i.e., 
Ft = {x±,X2, ■ ■ ■ ,Xp}, then it is possible to synchronize 
the entire lattice on that orbit. 

As an application, we shall consider a logarithmic map 
[TTj ] f(x) = b — In | x | , x G (—00,00), as local chaotic 
dynamics in both Eq. (JTJ and Eq. J2J ■ This map exhibits 
robust chaos, with no periodic windows and no separated 
chaotic bands, in the parameter interval b G [— 1, 1], 

In order to characterize the collective states in both 
systems, we calculate the mean field 



1 N 

<=1 



(7) 



Figures l(a)-l(c) show bifurcation diagrams of h t as a 
function of the coupling drive parameter e 2 for lattices 
driven with different forms of F t . For each value of e 2 , 
ht was calculated at each time step during a run start- 
ing from random initial conditions on the local maps, 
uniformly distributed on the interval [—8,4], after dis- 
carding the transients. Figure 1(a) shows the bifurcation 
diagram of h t vs. e 2 for a lattice driven with a constant 
term F t — x\ = —0.855, where x\ = f(x\) is the unsta- 
ble fixed point of the local logarithmic map for b = —0.7. 
In the region labeled by PS (periodic synchronization), 
ht becomes equal to x±, indicating that the elements in 
the lattice are synchronized at this stationary value. The 



range of e 2 for which synchronization is observed corre- 
sponds to the range predicted by the stability condition 
Eq. ©, for the unstable fixed point of / with p = 1 . 
Outside the region of stationary synchronization, other 
types of collective behaviors can be observed in the bi- 
furcation diagram of Fig. 1(a). In the left region labeled 
by CPB (collective periodic behavior), the mean field of 
the system driven with constant F t = x\ experiences a 
collective period-2 motion, although the local elements 
are chaotic and desynchronized. A collective fixed point 
occurs on the right CPB region. There is also a region, 
labeled by T (turbulence) , where h t follows a normal sta- 
tistical behavior around a mean value with fluctuations 
reflecting the averaging of N completely desynchronized 
chaotic elements [l2j. Figure 1(b) shows h t vs. e 2 for 
the lattice driven periodically with F t = {ii,a; 2 }, where 
xt = 0.18, 5 2 = —2.44 are the points of the unstable 
period-2 orbit of the local logarithmic map for b = —0.7. 
In the region labeled by PS, the mean field coincides with 
this unstable periodic orbit, indicating that the lattice is 
synchronized on this orbit. This range of e 2 is predicted 
by condition Eq. (jfJJl for the unstable periodic orbit of 
/ with p = 2 for the parameter value b = —0.7. Peri- 
odic collective behaviors with period 2 and period 4 in 
ht arise in the regions labeled by CPB. Figure 1(c) cor- 
responds to a chaotic driving, with F t = — 0.7 + In |art|. 
Chaotic synchronization of the system occurs in the re- 
gion labeled by CS and it is also predicted by Eq. |JBJ). In 
this CS region, h t = F t = -0.7 + In |jc t | = x\, Mi. After 
crossing the boundaries of the CS region, the collective 
states described by h t take the form of chaotic bands. 
These states are labeled CCB (collective chaotic bands) 
and they consist of the motion of chaotic elements that 
maintain some coherence. 

Synchronized states, x\ — Xt, Vz, can also emerge in 
the autonomous coupled map system Eq. I n order 
to study the stability of these states, we express Eq. 
in vector form as 



x f+ i 



(l- e2 )/+|L+|G 



f(xt), 



(8) 



where the local connectivity matrix L was defined in 
Eq. J2J) and G is the N x N global connectivity matrix 
having all its components equal to 1. 

The linear stability analysis of synchronized states in 
the autonomous system yields 



<1, j = 0,l,-,N-l, 



(9) 

where Vj are the eigenvalues of the local coupling ma- 
trix L defined above, and jj — 5 jN correspond to the 
set of eigenvalues of G, with the zero eigenvalue having 
(JV — l)-fold degeneracy. The range of e 2 for which syn- 
chronization takes place is 



l-2 ei sin 2 (^; 



A <e 2 < l-2e lS in 2 (^- 



(10) 



3 



T 


CPB 


PS 


' (a) 

CPB 


< 


I 






( 







-0.2 0.2 0.4 0.6 0.8 1 



ht 



— iii CPB /* 


] PS 


1 (b) 

CPB 


- 




/ 










\ 

\ 




1 


U 





0.2 


0.4 0.6 

£2 


0.8 




• CCB ,,. 






: If 


■ ■ ■ i CCB; ■ 












'. 1 ' ■ . . 


).2 





0.2 


0.4 0.6 

£2 


0.8 



FIG. 1: Bifurcation diagrams of ht vs. 62 for the driven lat- 
tice, Eq. JTJ, with TV = 10 4 , t\ = 0.5, 6 = -0.7. Labels 
indicate PS (periodic synchronization), CPB (collective peri- 
odic behavior), CS (chaotic synchronization), CCB (collective 
chaotic bands), and T (turbulence), (a) Ft = {xi = —0.855}. 
(b) Ft = {xi = 0.18, x 2 = -2.44}. (c) F t = -0.7 + In \x t \. 



which is the same condition for stability of synchronized 
states in the driven lattice, Eq. However, the unsta- 
ble periodic orbits of the local map / can not be synchro- 
nized in the autonomous system because they correspond 
to unstable synchronized states in this system. 

The eigenvector corresponding to j = is homoge- 
neous for both matrices L and G. Thus perturbations of 
Xj along this eigenvector do not destroy the coherence, 
and the stability condition associated with j = is ir- 
relevant for a synchronized state. The other eigenvalues 
corresponding to j ^ are not homogeneous. Thus, con- 



dition Eq. ) with j ^ define regions in the space of 
parameters where all the stable synchronized states can 
be observed. 

Figure 2 shows the boundary curves given by Eq. on 
the parameter plane (£2, ei), in the limit of large N. The 
label CS indicates where a chaotic synchronized state is 
stable in the autonomous system. After crossing the sta- 
bility boundaries of the CS state, the autonomous system 
exhibits spatiotemporal patterns corresponding to short 
wave (SW) and long wave (LW) modes, as indicated in 
Fig. 2. Note that the stability condition Eq. © does not 
hold for £2 = and large N; i.e., chaotic synchronization 
for large system size N can not take place in coupled 
map lattices possessing only local couplings, which is a 
well known fact 01 . Thus, the addition of a global in- 
teraction allows the emergence of chaotic synchronization 
in a large locally coupled system. Similarly, an external 
uniform driving can be used to induce chaotic synchro- 
nization in a locally coupled map network. 
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FIG. 2: Stability boundaries of the chaotic synchronized state 
CS of the autonomous system Eq. on the parameter plane 
(e 2 ,ei), with N = 10 4 , b = -0.7. Re gions corresponding to 
short wave (SW) and long wave (LW) patterns emerging from 
the CS state are indicated. 

Figure 3(a) shows the bifurcation diagram of ht for 
the autonomous system, Eq. J2J), as a function of the 
coupling parameter £2- In contrast to the behavior dis- 
played by the driven lattice, no regions of synchroniza- 
tion of unstable periodic orbits of the local dynamics in 
the space of parameters are observed in the autonomous 
system, as expected. The range of chaotic synchroniza- 
tion (CS) corresponds to the same range of £2 for chaotic 
synchronization in Fig. 1(c) for the driven lattice. Note, 
however, that beyond the region CS in Fig. 3(a), there 
are collective states emerging in the autonomous system, 
such as CPB (collective periodic behavior), and T (turbu- 
lence), that do not appear in the corresponding diagram 
of the driven lattice, Fig. 1(c). By varying the coupling 
strengths, various spatial patterns can be realized in the 
autonomous system. These patterns correspond to short 
and to long wave modes and can be regarded as general- 
ized Turing patterns that emerge when a global synchro- 
nized state becomes unstable [l4(- Figure 3(b) shows the 
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FIG. 3: Bifurcation diagrams of h t for the autonomous sys- 
tem, for b = -0.7, TV = 10 4 . (a) ht vs. e 2 , with fixed ei = 0.5. 
(b) /it vs. ei, with fixed £2 = 0.5. Labels are defined in Fig. 1 



bifurcation diagram of h t vs. e\, with e 2 fixed. 

The collective states observed in the bifurcation dia- 
grams of Fig. 3, as well as those in Fig. 1, have been 
checked for system-size effects. When the lattice size N 
is increased, the segments in the CPB regions shrink, but 
the amplitudes of the collective periodic motions do not 
decrease, indicating that the global periodic attractors 
become better defined in the large system limit. As a 
consequence, the variance of the fluctuations of h t itself 



do not decrease as TV -1 with increasing N, but rather 
it saturates at some constant value related to the ampli- 
tude of the collective periodic motion. These states of 
ordered evolution in macroscopic quantities such as ht, 
arising from local chaos in spatiotemporal_svstems, have 
been called nontrivial collective behavior [l5j ■ The states 
of collective chaotic bands (CCB) are also manifestations 
of nontrivial collective behavior, since the variance of h t 
for those states neither follows a regular statistical be- 
havior. In contrast, the variance of the mean field in the 
turbulent states (T) appears to decrease as iV -1 with 
increasing N, obeying normal statistical behavior. 

In summary, we have studied the analogy between the 
chaotic synchronized states emerging in forced spatiotem- 
poral systems and in autonomous dynamical systems pos- 
sessing global interactions in the context of coupled map 
lattices. By using a model of a one-dimensional coupled 
map lattice subjected to a uniform external drive, we 
have shown that both, synchronization of unstable peri- 
odic orbits of the local maps and chaotic synchronized 
states, can be induced in the driven lattice. The external 
drive acts as a control mechanism for stabilizing unsta- 
ble periodic orbits of the local maps. We showed that 
the synchronization behavior of the driven lattice can be 
compared, under some circumstances, with that of an au- 
tonomous coupled map system possessing a similar local 
coupling plus an additional global interaction that acts as 
a global feedback. It was found that the chaotic synchro- 
nized states in both systems are analogous; however, the 
autonomous system does not exhibit synchronization of 
unstable local periodic orbits. The collective states aris- 
ing when the chaotic synchronized state becomes unsta- 
ble are, in general, different in these systems. Although 
we have considered one-dimensional diffusive couplings, 
expressed by the matrix L, the analogy between a uni- 
form external drive and a global interaction can also be 
applied to other networks of coupled maps whose con- 
nectivity may be represented in matrix form. 
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